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Abstract 

We analyze the polarization observables of J ftp photoproduction at next-to- 
leading order (NLO) within the factorization formalism of nonrelativistic quantum 
chromodynamics (NRQCD). This is the first NLO study of heavy-quarkonium polar- 
ization including the full relativistic corrections due to the intermediate Vjq 8 ' , 3 S?, 
and 3 Pj 8 ^ color-octet (CO) states in the worldwide endeavor to test NRQCD fac- 
torization at the quantum level. We present theoretical predictions in the helicity, 
target, and Collins-Soper frames of DESY HERA, evaluated using the CO long- 
distance matrix elements previously extracted through a global fit to experimental 
data of unpolarized J ftp production, and confront them with recent measurements 
by the HI and ZEUS Collaborations. We find the overall agreement to be satisfac- 
tory, but the case for NRQCD to be not as strong as for the J/ip yield. 

PACS numbers: 12.38.Bx, 13.60.Le, 13.88.+e, 14.40.Pq 



The test of NRQCD factorization [T] has been identified to be among the most exigent 
milestones on the roadmap of heavy-quarkonium physics at the present time [2] . Quarko- 
nia are systems consisting of a quark and its antiparticle bound by the strong force, among 
which charmonium cc and bottomonium bb are considered heavy. The J ftp meson, the 
lowest-lying cc state of spin one, which was simultaneously discovered at the Brookhaven 
National Laboratory [3] and the Stanford Linear Accelerator Center [I] in November 1974 
(The Nobel Prize in Physics 1976), provides a particularly useful laboratory for such a 
test because it is copiously produced at all high-energy particle colliders, owing to its 
relatively low mass, and particularly easy to detect experimentally. In fact, sharing the 
total-angular-momentum, parity, and charge-conjugation quantum numbers J PC = 1 
with the photon, it can decay to e + e~ and pairs producing spectacular signatures 

in the detectors, the branching fraction of either decay channel being as large as about 
6% 0. 

In fact, the NRQCD factorization formalism pQ provides a rigorous theoretical frame- 
work for the description of heavy-quarkonium production and decay. This implies a sepa- 
ration of process-dependent short-distance coefficients, to be calculated perturbatively as 
expansions in the strong-coupling constant a s , from supposedly universal long-distance 
matrix elements (LDMEs), to be extracted from experiment. The relative importance of 
the latter can be estimated by means of velocity scaling rules, which predict each of the 
LDMEs to scale with a definite power of the heavy-quark (Q = c, b) velocity v in the limit 
!)<1. In this way, the theoretical predictions are organized as double expansions in a s 
and v. A crucial feature of this formalism is that the QQ pair can at short distances be 
produced in any Fock state n = 2S+1 l}-f with definite spin S, orbital angular momentum 
L, total angular momentum J, and color multiplicity a = 1,8. In particular, this for- 
malism predicts the existence of intermediate CO states in nature, which subsequently 
evolve into physical, color-singlet (CS) quarkonia by the nonperturbative emission of soft 
gluons. In the limit v — > 0, the traditional CS model (CSM) is recovered in the case of 
S- wave quarkonia. In the case of J/ip production, the CSM prediction is based just on 
the 3 Si CS state, while the leading relativistic corrections, of relative order 0(v 4 ), are 
built up by the 1 S [ o\ 3 S?, and 3 P| 8] (J = 0, 1, 2) CO states. 

The CSM is not a complete theory, as may be understood by noticing that the NLO 
treatment of P-wave quarkonia is plagued by uncanceled infrared singularities, which are, 
however, properly removed in NRQCD. This conceptual problem cannot be cured from 
within the CSM, neither by proceeding to higher orders nor by invoking kx factorization 
etc.. As it were, NRQCD factorization, appropriately improved at large transverse mo- 
menta Pt by systematic expansion in powers of tuq/p^ [6], is the only game in town, 
which makes its experimental verification such a matter of paramount importance and 
general interest [2]. 

The present status of testing NRQCD factorization in charmonium production is as 
follows. Very recently, NRQCD factorization has been consolidated at NLO by a global fit 
[7] to all available high-quality data of inclusive unpolarized J ftp production, comprising 
a total of 194 data points from 26 data sets collected by 10 experiments at 6 colliders, 
namely by Belle at KEKB; DELPHI at LEP II; HI and ZEUS at HERA I and II; PHENIX 
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at RHIC; CDF at Tevatron I and II; and ATLAS, CMS, ALICE, and LHCb at the LHC. 
This fit successfully pinned down the three CO LDMEs in compliance with the velocity 
scaling rules, establishing their universality, and yielded an overall description of the data 
well within the theoretical uncertainties; appreciable deviations arose only in the case of 
two-photon scattering, where the useable data comprises only 16 events and has not been 
confirmed by any of the other three LEP II experiments, however. On the other hand, the 
NLO CS predictions were found to significantly undershoot all the measurements, except 
for the single data point of e + e~ annihilation. 

In contrast to the J /if) yield, NRQCD interpretations of J/if> polarization measure- 
ments have so far been exhibiting a rather confusing pattern [2] , presumably because the 
theoretical status is much less advanced there. In fact, complete NRQCD predictions for 
J /if) polarization observables so far only exist at LO. At NLO, the CSM predictions for 
direct photoproduction [8f9] and hadroproduction [10| as well as the 1 Sq and 3 Sj 8 ' contri- 
butions to hadroproduction [II] . which may be obtained using standard techniques, are 
known. The NLO calculation of 3 Pf^ contributions, which are expected to be significant, 
is far more intricate because the applications of the respective projection operators to the 
short- distance scattering amplitudes produce particularly lengthy expressions involving 
complicated tensor loop integrals and exhibiting an entangled pattern of infrared singu- 
larities. This technical bottleneck is overcome here for the first time for J /if) polarization 
observables. 

Recent high-quality measurements by the HI [12] and ZEUS [13] Collaborations at 
HERA provide a strong motivation for us to start by studying photoproduction, where 
the incoming leptons interact with the protons via quasi-real photons, of low virtuality 
Q 2 = —p 2 , and are deflected under small angles. Such quasi- real photons participate 
in the hard scattering either directly or via partons into which they fluctuate (resolve) 
intermittently. However, resolved photoproduction is greatly suppressed, to the level of 
1% [7], by the cut z > 0.3 (0.4) applied by HI [12] (ZEUS [13]) and is thus neglected 
here. Here, z = {pj/^> ■ P P )/(p-y • P P ), with pj/^, p y , and p p being the four-momenta 
of the J /if) meson, photon, and proton, respectively, denotes the inelasticity variable, 
which measures the fraction of photon energy transferred to the J / if) meson in the proton 
rest frame. Another important variable of photoproduction is the jp invariant mass, 



The polarization of the J/if> meson is conveniently analyzed experimentally by mea- 
suring the angular distribution of its leptonic decays, which is customarily parametrized 
using the three polarization observables A, fi, and u, as [H] 



where 9 and are respectively the polar the azimuthal angles of / + in the J /if) rest frame. 
This definition, of course, depends on the choice of coordinate frame. Among the most 
frequently employed coordinate frames are the helicity (recoil), Collins-Soper, and target 



W = ^( P ,+P P ) 2 . 



dT(J/iP -> /+/-) 

d cos 9 d<f> 



oc 1 + A cos 2 9 + n sin(20) cos cf> 
+ ^sin 2 6>cos(20), 



(1) 
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(O-W(isf)) (3.04 ± 0.35) x 10- 2 GeV 3 
lo J ^( 3 Sf ] )) (1.68 ± 0.46) x 10- 3 GeV 3 
(Qj/i>(3pW)) (-9.08 ± 1.61) x IP' 3 GeV 5 



Table 1: J/tp NLO CO LDMEs corrected for feed-down [7J. 

frames, in which the polar axes point in the direction of — (p y + p p ), Py/\Py \ — Pp/\Pp\, an d 
— p p , respectively. The values A = 0, +1, — 1 correspond to unpolarized, fully transversely 
polarized, and fully longitudinally polarized J/ij) mesons, respectively. 
On the theoretical side, we have 

_ dan - dapo 

dau + daoo' 

V^Re da±o 

dan + daoo ' 
2d(7i-i 

v • = 3 j — , (2) 

dan + da 00 

where da^j, with i,j = 0, ±1 denoting the z component of S, is the ij component of the 
ep — > J/ip + X differential cross section in the spin density matrix formalism. Invoking the 
Weizsacker- Williams approximation and the factorization theorems of the QCD parton 
model and NRQCD [TJ, we have 

<H = £ / dxdyf, /e (x)f k/p (y)(0 J ^[n]) 

k,n 

x ^-dPS pijijk -> cc[n] + X), (3) 

where / 7 / e (x) is the photon flux function, fk/ P (y) the parton distribution function (PDF) 
of parton k = g } q,q with q = u, d,s, {0 J M[n]) are the LDMEs, s = (p^+Pk) 2 , and dPS 
is the phase space measure of the outgoing particles. The spin density matrix elements 
of the partonic cross sections, Pij(jk —> cc[n) + X), are averaged (summed) over the 
spins and colors of the incoming (outgoing) particles, keeping i and j fixed for the cc 
pair in Fock state n. The quantities pij are evaluated by applying polarization and color 
projectors similar to those listed in Ref. (15] to the squared QCD matrix elements of open 
cc production. For n = 3 Si~\ 3 Sf\ 3 Pf\ the cc[n] spin polarization vectors e(i) appearing 
in p^ are replaced by their explicit expressions [16] . In the case of n = 3 Pf \ for which 
S = L = 1, the z components of L are summed over. Forn = 1 4 8] ,p 11 and Poo are each 
set to one third of the squared matrix element, and pio = Pi-i = 0. For space limitation, 
we refrain from presenting here more technical details, but refer the interested reader to 
a forthcoming publication. 

We now describe the theoretical input for our numerical analysis. In all our NRQCD 
calculations, we use the CO LDME set extracted in Ref. [7] after subtracting from the 
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Figure 1: (color online) NLO NRQCD predictions (solid lines) for A and v as functions 
of px and z in the helicity, Collins-Soper, and target frames including theoretical uncer- 
tainties (shaded/yellow bands) compared to HI [12] and ZEUS [13] data. For compari- 
son, also the NLO CSM (dot-dashed lines) predictions including theoretical uncertainties 
(hatched/blue bands) as well as the LO NRQCD (dashed lines) and LO CSM (dotted 
lines) ones are shown. 



data fitted to the estimated contributions due to feed-down from heavier charmonia. For 
the reader's convenience, these values are listed in Table [TJ For consistency, we also 
adopt the residual input from Ref. [7]. In particular, we choose the CS LDME to be 
(O j ^( 3 S?)} = 1.32 GeV 3 [TF] and the charm-quark mass, which we renormalize according 
to the on-shell scheme, to be m c = 1.5 GeV, adopt the values of the electron mass m e and 
the electromagnetic coupling constant a from Ref. [5], and use the one-loop (two-loop) 
formula for oi n/ ^(/i r ), with nf = 4 active quark flavors, at LO (NLO). As for the proton 
PDFs, we use the CTEQ6L1 (CTEQ6M) set [18J at LO (NLO), which comes with an 
asymptotic scale parameter of Aq^. d = 215 MeV (326 MeV). We evaluate the photon 
flux function using Eq. (5) of Ref. [19J. Our default choices for the MS renormalization, 



5 



factorization, and NRQCD scales are fj, r = /if = ttit and //a = n^c, respectively, where 
rriT = \jp\ + 4m^ is the J /if) transverse mass. The bulk of the theoretical uncertainty is 
due to the lack of knowledge of corrections beyond NLO, which are estimated by varying 
fj, r , [if, and /ia by a factor 2 up and down relative to their default values. In our NLO 
NRQCD predictions, we must also include the errors in the CO LDMEs. To this end, 
we determine the maximum upward and downward shifts generated by independently 
varying their values according to Table [I] and add the resulting half-errors in quadrature 
to those due to scale variations. 

In Fig. [H we compare our NLO NRQCD predictions for A and v as functions of pt 
and z, evaluated from Eq. (j2J) with the respective differential cross sections inserted on 
the r.h.s., with the measurements by HI [12] in the helicity and Collins-Soper frames and 
by ZEUS [13] in the target frame. The HI data was taken during the years 2006 and 
2007, and corresponds to an integrated luminosity of 165 pb -1 , while the ZEUS analysis 
covers all data collected from 1996 through 2007, corresponding to 430 pb _1 . At HERA, 
27.5 GeV electrons or positrons were colliding with 820 GeV (920 GeV) protons before 
(since) 1998. As the admixture of 820 GeV protons in the ZEUS data sample is negligible, 
we take the cm. energy to be 318 GeV also there. We adopt the experimental acceptance 
cuts, indicated in each of the six frames of Fig. [fl except for the pr distribution by ZEUS 
in Fig. ED^c). Unfortunately, ZEUS did not impose any upper z cut, which poses two 
problems on the theoretical side. On the one hand, in the kinematic endpoint region, at 
z « 1, where the scattering becomes elastic, the cross section is overwhelmed by diffractive 
J /if) production, the treatment of which lies beyond the scope of our paper. On the other 
hand, the NRQCD expansion in v breaks down in the limit z — > 1, so that our fixed-order 
calculation becomes invalid. We avoid these problems by introducing the cut z < 0.95, 
accepting that the comparison with the ZEUS data then has to be taken with a grain of 
salt. 

For comparison, also the LO NRQCD as well as the LO and NLO CSM predictions 
are shown in Fig. [TJ In order to visualize the size of the NLO corrections to the hard- 
scattering cross sections, the LO predictions are evaluated with the same LDMEs. We 
observe that, in all the cases considered, the inclusion of the NLO corrections has a 
considerably less dramatic effect in NRQCD than in the CSM, where the normalizations 
and shapes of the various distributions are radically modified. This indicates that the 
perturbative expansion in a s converges more rapidly in NRQCD than in the CSM. Looking 
at the \(pt) distributions in Figs. [Tfa)-(c), we notice that NRQCD predicts large-pr J/if> 
mesons to be approximately unpolarized, both at LO and NLO, which is nicely confirmed 
by the HI measurements in Figs. IU(a) and (b). However, the ZEUS measurement in 
Fig. H](c), which reaches all the way up to z = 1, exhibits a conspicuous tendency towards 
transverse polarization, which might well reflect the notion that diffractively produced 
vector mesons prefer to be strongly transversely polarized at z ~ 1 [20]. Comparing the 
NLO NRQCD and CSM predictions in the three different frames, we conclude that the 
Colins-Soper frame possesses the most discriminating power. As expected, the theoretical 
uncertainties, which are chiefly due to scale variations, steadily decrease as the value of 
Pt increases, which just reflects asymptotic freedom. By the same token, the theoretical 



6 



uncertainties in the z distributions in Figs.[Tfd)-(e), which are dominated by contributions 
from the pt region close to the lower cut-off at Pt = 1 GeV, are quite sizable, which makes 
a useful interpretation of the experimental data more difficult. 

At this point, we compare our results with the theoretical literature. We agree with 
the LO NRQCD formulas for Pij(jk — > cc[n] + k) listed in Appendix B of Ref. [16J. We 
are able to nicely reproduce the NLO CSM results for A and v as functions of pt and z 
shown in Fig. 2 of Ref. [8] and Fig. 2 of Ref. [9] if we adopt the theoretical inputs specified 
there. The differences between those NLO CSM results and the respective results in our 
Fig. [T] are due to the use of different theoretical inputs. A similar statement applies to the 
LO NRQCD results graphically displayed in Refs. [T2"fl3|ll6] . which are evaluated using 
CO LDMEs obtained from LO fits to Tevatron I data. 

In contrast to the unpolarized J / ip yield, where the most precise world data uniformly 
and vigorously support NRQCD and distinctly disfavor the CSM at NLO [7J, the situation 
seems to be less obvious for the J/ip polarization in photoproduction, as a superficial 
glance at Fig. [T] suggests. However, detailed investigation reveals that the overall x 2 value 
of all the HI and ZEUS data in Fig. [TJw.r.t. the default NLO predictions is reduced by 
more than 50% as the CO contributions are included, marking a general trend towards 
continued verification of NRQCD factorization. Unfortunately, this is where the legacy of 
HERA, which was shut down in 2007, ends. With the help of the proposed lepton-proton 
collider LHeC at CERN, polarized J/ip photoproduction could be studied more precisely 
and up to much larger values of px- Fortunately, measurements of J/ip polarization have 
also been performed in hadroproduction at the Tevatron and will be carried on at the 
LHC for many years. This is arguably the last frontier in the international endeavor to 
test NRQCD factorization in charmonium physics. 

This work was supported in part by BMBF Grant No. 05H09GUE and HGF Grant 
No. HA 101. 
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